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Abstract 

I present the highlights of a recent study of the effective QCD phase diagram on the temperature T and quark chemical potential 
fi plane, where the strong interactions are modeled using the linear sigma model coupled to quarks. The phase transition line is 
found from the effective potential at finite T and ji taking into account the plasma screening effects. We find the location of the 
critical end point (CEP) to be jT^) ~ (1.2,0.8), where is the (pseudo)critical temperature for the crossover phase 

transition at vanishing /i. This location lies within the region found by lattice inspired calculations. Since the linear sigma model 
does not exhibit confinement, I argue that the location is due to the proper treatment of the plasma screening effects and not to the 
size of the confining scale. I also comment on the extension of this study to determine the dependence of the CEP’s location on the 
strength of an external magnetic field. 
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1. Introduction - 

The different phases in which matter made up of 
quarks and gluons arranges itself depends on the tem¬ 
perature and density, or equivalently, on the temperature 
and chemical potentials. The representation of the QCD 
phase diagram is thus two dimensional. This is custom¬ 
ary plotted with the light-quark chemical potential /r as 
the horizontal variable and the temperature T as the ver¬ 
tical one. ju is related to the baryon chemical potential 
fiB by jUb = 3fi. 

Most of our knowledge of the phase diagram is re¬ 
stricted to the ju = 0 axis. The phase diagram is, by and 
large, unknown. For physical quark masses and ju = 0, 
lattice calculations have shown 111] that the change from 
the low temperature phase, where the degrees of free¬ 
dom are hadrons, to the high temperature phase de¬ 
scribed by quarks and gluons, is an analytic crossover. 
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Th e phas e transition has a dual natur e : on th e on e hand 

the color-singlet hadrons break up leading to deconfined 
quarks and gluons; this is dubbed as the deconfinement 
phase transition. On the other hand, the dynamically 
generated component of quark masses within hadrons 
vanishes; this is referred to as chiral symmetry restora¬ 
tion. 

Lattice calculations have provided values for the 
crossover (pseudo)critical temperature Tc for p - 0 
and 2 - 1-1 quark flavors using different types of im¬ 
proved rooted staggered fermions. The MILC collab¬ 
oration obtained Tc = 169(12)(4) MeV. The RBC- 
Bielefeld collaboration reported Tc - 192(7)(4) MeV. 
The Wuppertal-Budapest collaboration has consistently 
obtained smaller values, the latest being Tc = 147(2)(3) 
MeV. The HotQCD collaboration has computed Tc - 
154(9) MeV and more recently Tc = 155(1)(8) 

MeV H. The differences could perhaps be attributed 
to different lattice spacings. 

Although the above picture presented by lattice QCD 
cannot be easily extended to the case p 0 due to 
the sign problem, some mathematical extensions of lat¬ 
tice techniques |^ as well as Schwinger-Dyson equa- 
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tions 101, can be employed to explore all the phase dia¬ 
gram. 

A number of different model approaches indicate that 
the transition along the fi axis, at T = 0, is strongly first 
order ||3l. Since the first order line originating at T = 0 
cannot end at the /i = 0 axis which corresponds to the 
starting point of the cross-over line, it must terminate 
somewhere in the middle of the phase diagram. This 
point is generally referred to as the critical end point 
(CEP). The mathematical extensions of lattice tech¬ 
niques place the CEP in the region (//“’’/Tc, r™’’/7’c) ~ 
(1.0- 1.4,0.9-0.95) & 

In the first of Refs. IJ], it is argued that the theo¬ 
retical location of the CEP depends on the size of the 
confining length scale used to describe strongly inter¬ 
acting matter at finite density/temperature. This argu¬ 
ment is supported by the observation that the models 
which do not account for this scale 100 produce ei¬ 
ther a CEP closer to the axis (/r™’’/7’c and IT c 

larger and smaller, respectively) or a lower Tc 10 than 
the lattice based approaches or the ones which consider 
a finite confining length scale. Given the dual nature 
of the QCD phase transition, it is interesting to explore 
whether there are other features in models which have 
access only to the chiral symmetry restoration facet of 
QCD that, when properly accounted for, produce the 
CEP’s location more in line with lattice inspired results. 

An important clue is provided by the behavior of the 
critical temperature as a function of an applied magnetic 
field. Lattice calculations have found that this temper¬ 
ature decreases when the field strength increases il2l - 
0 . It has been recently shown that this phenomenon, 
dubbed inverse magnetic catalysis, can be obtained in 
models, such as the Abelian Higgs model or the lin¬ 
ear sigma model with quarks, which show only chiral 
symmetry restoration and lack confinement. This re¬ 
sult is a consequence of the decrease of the coupling 
constants with increasing field strength. The novel fea¬ 
ture implemented in these calculations is the handling 
of the screening properties of the plasma, which effec¬ 
tively makes the treatment go beyond the mean field ap¬ 
proximation 00 and allows to consider the ther- 
momagnetic modifications of the coupling constants at 
lowest order within the same calculation. Screening is 
also important to obtain a decrease of the coupling con¬ 
stant with the magnetic field strength in QCD in the 
Hard Thermal Loop approximation 11711 . 

It therefore seems that properly accounting for the 
plasma screening effects in effective models allows to 
obtain both a CEP’s location in line with lattice in¬ 
spired techniques as well as inverse magnetic cataly¬ 
sis. A pertinent question is what happens to the CEP’s 


location once a magnetic field dependence is included 
in the analysis for p 4 0 and wether the nature of the 
phase transition for - 0 changes as the magnetic field 
strength increases. Recent lattice QCD calculations IjJ | 
show that at very high values of the magnetic field 
strength, inverse magnetic catalysis prevails and that the 
phase transition becomes first order at asymptotically 
large values of the magnetic field for /i = 0 (see also 
Ref. ( 01 ). In this work we explore the consequences 
of the proper handling of the plasma screening proper¬ 
ties in the description of the effective QCD phase dia¬ 
gram within the linear sigma model with quarks. We ar¬ 
gue that it is the adequate description of these properties 
which determines the CEP’s location. We find that for 
certain values of the model parameters, obtained from 
physical constraints, the CEP’s location agrees with lat¬ 
tice inspired calculations. We also give a preview of 
work in progress iItoI] that shows that when including 
the effects of a magnetic field in the calculation of both 
the effective potential as well as on the thermomagnetic 
dependence of the coupling constants, the CEP’s loca¬ 
tion moves toward smaller values of the chemical po¬ 
tential and lower temperatures and that above a certain 
value of the field strength the CEP reaches the T -axis 
and the phase transitions become first order, also in line 
with recent lattice results 11811 . Details of the calcula¬ 
tion that forms the basis of this work can be found in 

Ref. 0 . 


2. The linear sigma model with quarks 

We start from the linear sigma model coupled to 
quarks. It is given by the Lagrangian density 

- ^(cr^ + + i^yO^i// 

- g^{a- + ij5f-n)ip, ( 1 ) 

where ijj is an SU(2) isospin doublet, jf = (7Ti,n2,tto,) 
is an isospin triplet and cr is an isospin singlet. The 
neutral pion is taken as the third component of the pion 
isovector, and the charged pions as 7t± - (tti + 

in2)l2. The squared mass parameter a^ and the self¬ 
coupling A and g are taken to be positive. 

To allow for the spontaneous breaking of symmetry, 
we let the cr field develop a vacuum expectation value v 

cr — > cr H- V, (2) 

which can later be taken as the order parameter of the 
theory. After this shift, the Lagrangian density can be 
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rewritten as 
£ = 


■^crd^d^o--^ (3/lv^ - a^) < 


- —jidud'^n - (Av^ — a^) if + — 

- + *iA/ d^iff - gvif/tf/ + X; + X{, (3) 

where Xt and Xf are given by 


fb - _ 


(cr2 + {nyf 

+ An^n~{cr^ + + ;7r^;7r“) 


jy 


■f _ 


-gip(cr + iy^f ■ n)tp. 


(4) 


and describe the interactions among the fields cr, n and 
if/, after symmetry breaking. From Eq. ([3]) we see that 
the cr, the three pions and the quarks have masses 


'«/ 


3Av^ - a^, 
Av^ - a^, 
gv, 


(5) 


respectively. 

The one-loop effective potential for the linear sigma 
model with quarks including the plasma screening prop¬ 
erties encoded in the ring diagrams contribution has 
been calculated in detail for zero chemical potential in 
Refs. 1 3 23]. Such analyses show that inclusion of the 
ring diagrams renders the effective potential stable. 

When the fi is non-vanishing, the calculation of the 
effective potential is more complicated. Though the bo¬ 
son contribution remains the same, the fermion contri¬ 
bution has to be modified due to the chemical potential. 
The modification enters the calculation in two ways; in¬ 
directly into the boson self-energy and directly from its 
contribution to the effective potential. 

To one-loop order the fermion contribution to the 
effective potential in the imaginary time formalism of 
thermal field theory is given by 122 j 


Vf 


-- r 

ln(l 


d^k 


+ e 


(2;r)3 


[/3w -H ln(l 


+ e 




( 6 ) 


where /? = and cj - {!? + and the sum 

over the fermion Matsubara frequencies has been per¬ 
formed. The first term in Eq. (|6]l corresponds to the 
vacuum contribution whereas the second and third ones 
are the matter contributions. Note that the matter con¬ 
tribution is made out of separate quark and antiquark 
pieces due to the finite chemical potential. The vacuum 


contribution is well-known Il22l] and can be expressed, 
after mass renormalization as a function of the renor¬ 
malization scale yu. For the evaluation of the medium’s 
contribution in Eq. (|6ll we adapt the technique from 
Ref. l24|] to the present case. The main idea is to pro¬ 
duce a second-order differential equation in y^, where 
y - niflT, valid at high temperature with nif as the 
smallest of all scales, for the finite temperature part of 
the potential, which we denote by Vf, given in Eq. (|6ll 
with appropriate boundary conditions at y = 0, where 
the integrals can be analytically evaluated. The expres¬ 
sion for the effective potential is obtained by integrating 
this differential equation and using the given boundary 
conditions. Combining the vacuum contribution after 
mass renormalization with the finite temperature part 
we finally have Ei 


Vf-- 


1 


16;r^ 


In 


(47rr)^ 

Iff- 


+ lA” - -H 
^ \ 2 IttT 




2nTj 


- 32r'‘[L4(-e'^^^)-i-L4(-e“^^^)]}. (7) 

It can also be shown that the boson self-energy If 
computed for a finite chemical potential and in the limit 
where the masses are small compared to T, is given by 

AT^ 

NfNcg^T^ 


n = 


,( 8 ) 


where Nf - 2 and = 3 are the number of light flavors 
and colors, respectively. 

Choosing the renormalization scale as /i = e the 
effective potential up to the ring diagrams contribution 
is then given by 


y 


(eff) _ 


Cl ^ A 4 
■TT'' + 

2 4 


Zl 


m 


|^64;r^ 


In 


(AnTY 

2f 


- 2yE + 1 


2't'^ 


90 




Nc 

16 f 


Z 

f=u,d 


mfT 

24 


(AnTY 

2f 


^ \2 2nT 


In 




H- 1 


2nT 
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+ 8 mjT^ 

- 32 [L4(-e^^^) + Li4(-e~>^/^)]]. (9) 


In the limit when // —> 0, Eq. (|9]l becomes the expres¬ 
sion found in Refs. 3- In the same limit, Eq. (|8]l 
reduces to the well known expression for the self-energy 
at high temperature iQ. Equation (|3 represents the 
effective potential computed beyond the mean held ap¬ 
proximation that accounts for the leading screening ef¬ 
fects at high temperature. 


3. The phase diagram 


In order to hnd the values of the parameters A, g and 
a appropriate for the description of the phase transition, 
we note that when considering the thermal effects the 
boson masses are modihed since they acquire a thermal 
component. Eor ju = 0 they become 


mliT) 

ml(T) 


2t2 


3Av^ 


Av^ - ■ 


AT^ NfNcg^T 


2 6 
AT^ NfNcg^T^ 


( 10 ) 


At the phase transition, the curvature of the effective 
potential vanishes for v = 0. Since the boson thermal 
masses are proportional to this curvature, these also van¬ 
ish at V = 0. Erom any of the Eqs. (Uni), we obtain a 
relation between the model parameters at Tc given by 


^ A NfN,g2 

a — 1 r "X -1-■ 

\2 6 


( 11 ) 


Eurthermore, we can fix the value of a by noting from 
Eqs. (|5]l that the vacuum boson masses satisfy 


a — 


. - 3m1 


( 12 ) 


Since in our scheme we consider two-flavor massless 
quarks in the chiral limit, we take Tc - 170 MeV 02-511 
which is slightly larger than Tc obtained in Nj - 2 + 
1 lattice simulations. Also, in order to allow for a 
crossover phase transition for ji - 0 (which in our de¬ 
scription corresponds to a second order transition) with 
g, A ~ 1 we need that g^ > A. Since the effective po¬ 
tential is written as an expansion in powers of a IT we 
need that this ratio is smaller than 1. Erom Eqs. 
and ([12]) the coupling constants are proportional to nia- 
which, from the above conditions, restricts the analy¬ 
sis to considering not too large values of nia-. Since the 
purpose of this work is not to pursue a precise determi¬ 
nation of the couplings but instead to call attention to 
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Figure 1. EflFective QCD phase diagram computed for A = 0.2 and 
g = 0.71 obtained by considering = 300 MeV. For small values 
of the phase transition is second order. The order of the transition 
changes to first order for larger values of fi. The CEP is located at 
^cepit^ T^epit^^ ~ (1.2,0.8). 


the fact that the proper treatment of screening effects al¬ 
lows the linear sigma model to provide solutions for the 
CEP, we consider small values for mo-- Given that cr is 
anyhow a broad resonance, in order to satisfy the above 
requirements let us take for definitiveness mo- = 300 
namely, close to the two-pion threshold. Therefore, the 
allowed values for the couplings A and g are restricted 
by 


A ^ NfNcg^ 

2 6 


0.77. 


(13) 


Equation ([T3T l provides a relation between A and g. 
A possible solution consistent with the above require¬ 
ments is given as an illustration by d = 0.2, g - 
0.71. The corresponding phase diagram thus obtained 
is shown in Eig.[T] 

Note that for small ji the phase transition is second 
order. In this case the (pseudo)critical temperature is 
determined from setting the second derivative of the ef¬ 
fective potential in Eq. (|9]) to zero at v = 0. When yu 
increases, the phase transition becomes first order. The 
critical temperature is now computed by looking for the 
temperature where a secondary minimum for v ^4 0 is 
degenerate with a minimum at v = 0. In both of these 
cases, from the detailed analysis, we locate the position 
of the CEP as {^f^|Tc, T^^jTc) ~ (1.2,0.8), which is in 
the same range as the CEP found from lattice inspired 
analyses Note also that the phase transition curve is 
essentially flat close to the T axis. 
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4. Conclusions 

In conclusion, we have shown that it is possible to 
obtain values for the couplings that allow to locate the 
CEP in the region found by mathematical extensions of 
lattice analyses. Since the linear sigma model does not 
have confinement we attribute this location to the ade¬ 
quate description of the plasma screening properties for 
the chiral symmetry breaking at finite temperature and 
density. Magnetic field effects can be included in the de¬ 
scription i2(]ll both into the effective potential and into 
the behavior of the couplings, at lowest order. These last 
corrections lead to a decreasing of the couplings with 
the field strength. This decrease can be understood in 
general terms since the magnetic field produces a di¬ 
mensional reduction whereby the virtual particles that 
make up the vacuum are effectively constrained to oc¬ 
cupy Landau levels and thus restrict its motion to planes. 
This produces that charged virtual particles lie closer to 
each other and thus, because of asymptotic freedom, re¬ 
duce the strength of the interaction. This happens no 
matter how weak the external field may be. We have 
found also that as the field strength increases, the CEP’s 
location moves to lower values of p and of T and that in 
fact there is a value for this field where the CEP reaches 
the T -axis where the first order phase transitions remain 
for larger values of the field. These findings will be re¬ 
ported elsewhere shortly. We believe this description 
will play an important role in determining the location 
of the CEP also in QCD with and without magnetic 
fields. 
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